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Abstract 

Random Phase Approximation (RPA) provides a very convenient tool to study the ensembles of weakly interacting 
waves, commonly called Wave Turbulence. In its traditional formulation, RPA assumes that phases of interacting 
waves are random quantities but it usually ignores randomness of their amplitudes. Recently, RPA was generalised 
in a way that takes into account the amplitude randomness and it was applied to study of the higher momenta and 
probability densities of wave amplitudes. However, to have a meaningful description of wave turbulence the RPA 
properties assumed for the initial fields must be proven to survive over the nonlinear evolution time, and such a proof 
is the main goal of the present paper. We derive an evolution equation for the full probability density function which 
contains the complete information about the joint statistics of all wave amplitudes and phases. We show that, for any 
initial statistics of the amplitudes, the phase factors remain statistically independent uniformly distributed variables. 
If in addition the initial amplitudes are also independent variables (but with arbitrary distributions) they will remain 
independent when considered in small sets which are much less than the total number of modes. However, if the size 
^ of a set is of order of the total number of modes then the joint probability density for this set is not factorisable 
^ ^1 into the product of one-mode probabilities. In the other words, the modes in such a set are involved in a "collective" 
(correlated) motion. We also study new type of correlators describing the phase statistics. 



Introduction 



J>W^ave Turbulence (WT) is a common name for the fields of dispersive waves which are engaged in stochastic weakly 
'kJionlinear interactions over a wide range of scales. Plentiful examples of WT are found in oceans, atmospheres, plasmas 
jj,nd Bose-Einstein condensates 121 B El El • Roughly, there have been three major approaches to derive the WT 
Cd;heory, one based on a diagrammatic approach |S1 ^1 E| , the second based on cumulant expansions [21 E IllL [7] and the 
third one, the random phase approximation (RPA) 0E1E1- 

The diagrammatic approach was developed in a field theoretical spirit based on the Wyld's technique ^ . This method 
introduces an artificial Gaussian forcing for which a zero limit is taken at the end of the derivation. It is usually said that 
the statistical properties of this force (Gaussianity) do not affect the statistical properties of the resulting WT state which 
will be determined by the nonlinear properties only. However, such independence of the WT state on the statistics of the 
"seed" forcing is not obvious because the limit of small nonlinearity is taken before the limit of small force, i.e. the force 
remains much greater than the nonlinearity. In particular, when the nonlinearity parameter is strictly zero, the Wyld 
technique gives a Gaussian steady state which is clearly an artefact of this method because for linear systems statistics 
of the wave amplitudes remain the same as in the initial condition and, therefore, can be arbitrary. The question if any 
nonlinearity, no matter how small, can break this dependence of the steady state on the initial conditions still has not 
been answered in the literature. Thus, the diagrammatic approach, although a very efficient way to build the perturbation 
expansion, needs to be expanded to include non-Gaussian "seed" force in order to see to what extent the results are not 
sensitive to the force statistics. However, some elements of the Wyld technique will be used in the present paper, not as 
a complete description but rather as an auxiliary aid in writing out complicated terms. 
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The cumulant expansion approach differs from the other methods by working directly with the continuous Fourier 
transforms corresponding to the infinite coordinate space without introducing a finite box as an intermediate step. The 
main idea here is that, although the Fourier transform is ill-defined for the wave fields corresponding to homogeneous 
turbulence, it is well defined for the cumulants provided that the correlations decay rapidly enough in the coordinate 
space. The cumulant method is very elegant for describing the spectra and the multiple-point moments the points in 
which are not "fused" (i.e. all different). However, some important statistical quantities involve fused moments and they 
are hard (if at all possible) to define without introducing a finite box as an intermediate step. For example, one of such 
objects, (lal**), is important because it describes intensity of fluctuations of the fc-space distribution of energy |ap, namely 
5E = - (|a|2)2 (see ^). Furthermore, non-decaying in the x-space correlations tend to naturally develop over 

the nonlinear time ^] and it is not clear what wave fields these correlations correspond to within the cumulant approach. 

RPA approach has been by far most popular technique due to its clear intuitive content. However, this approach has 
occasionally been downgraded to just a convenient way of interpreting the results of a more rigorous technique based on the 
cumulant expansions. It happened because RPA, being widely used by physicists, had not been formulated rigorously. In 
particular, it is typically assumed that the phases evolve much faster than amplitudes in the system of nonlinear dispersive 
waves and, therefore, the averaging may be made over the phases only "forgetting" that the amplitudes are statistical 
quantities too (see e.g. Q]). This statement become less obvious if one takes into account that we are talking not about 
the linear phases ujt but about the phases of the Fourier modes in the interaction representation. Thus, it has to be 
the nonlinear frequency correction that helps randomising the phases ^J). On the other hand, for three- wave systems 
(considered in this paper) the period associated with the nonlinear frequency correction is of the same order in small 
nonlinearity e as the nonlinear evolution time and, therefore, phase randomisation cannot occur faster that the nonlinear 
evolution of the amplitudes. One could hope that the situation is better for 4-wave systems (not considered here) because 
the nonlinear frequency correction is still ~ but the nonlinear evolution appears only in the order. However, in order 
to make the asymptotic analysis consistent, such correction has to be removed from the interaction-representation 
amplitudes and the remaining phase and amplitude evolutions are, again, at the same time scale (now 1/e'*). This picture 
is confirmed by the numerical simulations of the 4-wave systems ^1 120] which indicate that the nonlinear phase evolves 
at the same timescale as the amplitude. Thus, to proceed theoretically one has to start with phases which are already 
random (or almost random) and hope that this randomness is preserved over the nonlinear evolution time. In most of 
the previous literature such preservation was assumed but not proven. The goal of this paper will be to study the extent 
to which such an assumption is valid. 

Another goal of this paper is to make RPA formulation more consistent by taking into account that both phases and the 
amplitudes are random variables. Indeed, even if one starts with a wavefield which has random phases but deterministic 
amplitudes, as it is typically done in numerical simulations, the amplitudes will get randomised because the nonlinear 
term producing their evolution contains (random) phase factors. Preliminary steps were recently done in I2()| where 
we assumed that all the phases and the amplitudes in the initial wavefield are random variables independent of each 
other and that the phase factors are uniformly distributed on the unit circle on the complex plane. We kept the same 
acronym RPA but re-interpreted it as "Random Phases and Amplitudes", reflecting the fact that, first, the amplitudes 
are also random and, second, that it is not an "approximation" but rather an assumed property of the initial field. Such 
a generalised RPA was used in to study the evolution of the higher moments of the Fourier amplitudes and in [21] 
to study their "one- mode" PDF. In fact, this form of RPA is more general than the cumulant approach because it can 
handle fields with long correlation lengths which appear to be important for intermittency |2()j . 

Of course, for such an analysis to be trustworthy one should prove that the RPA properties hold over the nonlinear 
time and not just for the initial fields. Such mathematical validation of the RPA method will be in the focus of the 
present paper. To do this we will have to study the full joint PDF which involves the complete statistical information 
about the system, including the multi-mode correlations of both the amplitudes and the phase factors. We will derive 
an evolution equation for such PDF and we will show that it is identical to the equation obtained for the excitations in 
anharmonic crystals originally obtained by Peierls ^H] and later reproduced by Brout and Prigogine ^Hl and Zaslavski 
and Sagdeev |17| . All these works were restricted to considering a quite narrow class of interaction Hamiltonians arising 
from a potential energy, i.e. depending on the coordinates but not momenta. These class does not include a large number 
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of interesting WT systems, e.g. the capillary, internal, Rossby and Alfven waves. It is remarkable, therefore, that the 
Peierls equation turns out to be universal for the general class of three- wave systems, as it is shown in the present paper. 
Further, we use this equation to validate an "essential" RPA formulation, i.e. approximate RPA which holds only up to 
a certain order in nonlinearity and discreteness, but which is sufficient for the WT closure. This validation gives RPA 
technique a status of a rigorous approach which, due to the simplicity of its premises, is a winning tool for the future 
theory of non-Gaussianity of WT, its intermittency and interactions with coherent structures. 

In addition to the mathematical validation of RPA, we will also develop WT further by considering new statistically 
important quantities. For a long time, describing and predicting the energy spectra was the only concern in WT theory. 
Recently, we presented a description of the higher order statistics of the one-point Fourier correlators in terms of their 
moments and PDF's. They describe the k-space "noise" , i.e. the fluctuations of the mode energy about its mean value 
given by the energy spectrum. We also showed PDF's have a long algebraic tail which indicates presence of intermittency 
in WT fields. The present paper deals with phases, and we will therefore introduce and study some new correlators which 
will allow to describe the phase statistics directly. Such a description will compliment the mathematical validation of the 
RPA because it yields to a physical answer on how initially correlated phases can get de-correlated in the first place. 

2 Fields with Random Phases and AmpUtudes. 

Let us consider a wavefield a(x, in a periodic cube of with side L and let the Fourier transform of this field be ai{t) 
where index IgZ'^ marks the mode with wavenumber ki = 2ttI/L on the grid in the d-dimensional Fourier space. For 
simplicity let us assume that there is a maximum wavenumber k^ax (fixed e.g. by dissipation) so that no modes with 
wavenumbers greater than this maximum value can be excited. In this case, the total number of modes is N = (kmax/'J^L)'^. 
Correspondingly, index / will only take values in a finite box, I G Bn C Z'^ which is centred at and all sides of which 
are equal to kmax/i^L — N^^^ . To consider homogeneous turbulence, the large box limit N ^ oo will have to be taken. ^ 
Let us write the complex ai as a; = Ai^i where Ai is a real positive amplitude and i/'; is a phase factor which takes 
values on 5^, a unit circle centred at zero in the complex plane. Let us define the A^-mode joint PDF 7^(^) as the 
probability for the wave intensities Af to be in the range (s/, si + dsi) and for the phase factors ipi to be on the unit-circle 
segment between and ^; -I- d^i for all I S Bn- In terms of this PDF, taking the averages will involve integration over all 
the real positive si's and along all the complex unit circles of all i^/'s. 



where notation /{A^jip} means that / depends on all Af^s and all V'z's in the set {Af,tpi;lEBN} (similarly, {s,^} means 
{sijipfj S Bn}, etc). The full PDF that contains the complete statistical information about the wavefield a(x, i) in the 
infinite x-space can be understood as a large-box limit 



i.e. it is a functional acting on the continuous functions of the wavenumber, Sk and S,k- In the the large box limit there 
is a path-integral version of 



The full PDF defined above involves all N modes (for either finite N or in the N oo limit). By integrating out all the 
arguments except for chosen few, one can have reduced statistical distributions. For example, by integrating over all the 

^It is easily to extend the analysis to the infinite Fourier space, kmax = oo. In this case, the full joint PDF would still have to be defined 
as a A'^ ^ oo limit of an TV-mode PDF, but this limit would have to be taken in such a way that both kmax and the density of the Fourier 
modes tend to infinity simultaneously. 




(1) 




(2) 
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angles and over all but M amplitudes, we have an "M-mode" amplitude PDF, 

^..,......,.M = ( n / dsi XI I \dU\\v^''\s,i}, (3) 

which depends only on the M amplitudes marked by labels ji, j2, • ■ • ,jM^BN- 
2.1 Definition of an ideal RPA field 

Following the approach of |19[ I20| . we now define a "Random Phase and Amplitude" (RPA) field. We say that the field 
a is of RPA type if it possesses the following statistical properties: 

1. All amplitudes Ai and their phase factors ifii are independent random variables, i.e. their joint PDF is equal to the 
product of the one-mode PDF's corresponding to each individual amplitude and phase, 



leBi. 



2. The phase factors V'/ are uniformly distributed on the unit circle in the complex plane, i.e. for any mode I 

Note that RPA does not fix any shape of the amplitude PDF's and, therefore, can deal with strongly non-Gaussian 
wavefields. Such study of non-Gaussianity and intermittency of WT was presented in |19[ 1201 and will not be repeated 
here. However, we will study some new objects describing statistics of the phase. 

In ^1120] RPA was assumed to hold over the nonlinear time. The main goal of this paper is to find out whether it is 
true that the RPA property survives over the nonlinear time and to what extent. We will see that RPA fails to hold in its 
pure form as formulated above but it survives in the leading order so that the WT closure built using the RPA is valid. 
We will also see that independence of the the phase factors is quite straightforward, whereas the amplitude independence 
is subtle. Namely, M amplitudes are independent only up to a 0{M/N) correction. Based on this knowledge, and leaving 
justification for later on in this paper, we thus reformulate RPA in a weaker form which holds over the nonlinear time 
and which involves Af-mode PDF's with M ^ N rather than the full A^-mode PDF. 

2.2 Definition of an essentially RPA field 

We will say that the field a is of an "essentially RPA" type if: 

1. The phase factors are statistically independent and uniformly distributed variables up to O(e^) corrections, i.e. 

^^^^{^.a = (^^^"^'"^4 [1 + 0(6^)], (4) 



where 



is the A^-mode amplitude PDF. 



2. The amplitude variables are almost independent is a sense that for each M <^ N modes the M-mode amplitude 
PDF is equal to the product of the one-mode PDF's up to 0{M/N) and o(e^) corrections, 

Vn,n,...,M = ^"^4"^ • ■ • Pm [1 + OiM/N) + 0{e^)]. (6) 
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3 Weak nonlinearity and separation of time scales 

Consider weakly nonlinear dispersive waves in a periodic box. Here we consider quadratic nonlinearity and the linear 
dispersion relations u>k which allow three-wave interactions. Example of such systems include surface capillary waves [S] 
112) . Rossby waves |13) and internal waves in the ocean |14j . In Fourier space, we have the following Hamiltonian equations, 

oo 
m,n— 1 

where a; = a{ki) is the complex wave amplitude in the interaction representation, ki — 2tiI/ L is the wavevector, L is the 
box side length, = Wfc, — cjfe^^ — LOk^ , loi — Wfc, is the wave frequency, is an interaction coefficient and e is a 
formal small nonlinearity parameter. 

In order to filter out fast oscillations at the wave period, let us seek for the solution at time T such that 27r/tj ^ T <C 
l/we^. The second condition ensures that T is a lot less than the nonlinear evolution time. Now let us use a perturbation 
expansion in small e, 

aKT)=a(°)+eaW+62aP^ (8) 

Substituting this expansion in Q we get in the zeroth order a^^^\T) — ai(0), i.e. the zeroth order term is time independent. 
This corresponds to the fact that the interaction representation wave amplitudes are constant in the linear approximation. 
For simplicity, we will write af^ (0) = ai , understanding that a quantity is taken at T = if its time argument is not 
mentioned explicitly. The first order is given by 

oc 

af\T) = -i ^ (V;J„,a„a„A^„(5^„+„ + 2V;™a,„a„A[;:^(5,"„) , (9) 

m,n— 1 

where Aj„„ — /J" e'"">"*(it = (e*'^™"^ — Here we have taken into account that af\T) — ai and a^i}\o) — 0. 

Iterating one more time we get 



where we used a^^\Q) — and introduced E{x, y) = A{x — y)e^y*dt. 



(10) 



4 Evolution of the multi-mode PDF 

In this section we will apply the approach of |19[ 1201 to derive the evolution equation for the multi-mode PDF via 
introducing a generating functional, performing a weak-nonlinearity expansion and statistical averaging aided by a new 
graphical technique. We are going to demonstrate the phase independence property. This will also prepare us to answer 
the question of the next section: to what extend the amplitudes are going to remain statistically independent? 
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4.1 Generating functional. 



Introduction of generating functionals often simplifies statistical derivations but it can be defined differently to suit a 
particular technique. For our problem, the most useful form of the generating functional is 

^^"'HA,^} = -l^(n e^'^'i^n^ (11) 
^ ' leBN 

where {A, /i} = {A;, £ Bn} is a set of parameters, \i<E.TZ and ^i^Z. 

where {^} = {/i; £ Z;l £ Bn}- This expression can be verified by considering mean of a function f{A^,^} using the 
averaging rule ^ and expanding / in the angular harmonics Tp™; m £ Z (basis functions on the unit circle), 

/{A^^} = ^.9{m,^} n V^rS (13) 

{m} leBN 

where {m} = {mi £ Z;l £ Bn} are indices enumerating the angular harmonics. Substituting this into with PDF 
given by H12|l and taking into account that any nonzero power of will give zero after the integration over the unit 
circle, one can see that LHS=RHS, i.e. that (|f 2(1 is correct. Now we can easily represent ((1211 in terms of the generating 
functional, 



(14) 



where stands for inverse the Laplace transform with respect to all A; parameters and {/i} = {/i; £ Z;l £ Bn} are 
the angular harmonics indices. 

By definition, in RPA fields all variables Ai and tpi are statistically independent and ^Ai's are uniformly distributed on 
the unit circle. Such fields imply the following form of the generating functional 

zW{A,m}-^(^'"^{A} n Sif^i)^ (15) 

where 

^(A^,a){A} = ( II e^'^') = Z(^){A,A*}U=o (16) 
leBN 

is an TV-mode generating function for the amplitude statistics. Here, the Kronecker symbol S{fj,i) ensures independence 
of the PDF from the phase factors ipi. As a first step in validating the RPA property we will have to prove that the 
generating functional remains of form 1(15(1 up to 1/A^ and 0{e^) corrections over the nonlinear time provided it has this 
form at t = 0. 



4.2 Asymptotic expansion of the generating functional. 

Let us first obtain an asymptotic weak-nonlinearity expansion for the generating functional Z{A,/i} exploiting the sepa- 
ration of the linear and nonlinear time scales. ^ To do this, we have to calculate Z at the intermediate time t = T via 
substituting into it ajiT) from (jS)) and retaining the terms up to O(e^) only. This calculation is given in the Appendix 
and the result of it is: 

Z{\ /i, T } = X{\, ^i, T} + X{X, -M, T} (17) 

^Hereafter we omit superscript (A'^) in the Af-mode objects if it does not lead to a confusion. 



6 



with 

where 
Jl 

J-i 



ll'IKA^ 



2KT 



\ I 3 



^ _ 

2 '"4|a(")|4 V 2 



(18) 

(19) 
(20) 
(21) 
(22) 



9\ll'^i 2^^3^k(aj a- + (Aj + J— — -(a^ - a^^ ja^. a^. \ , (2^) 

^ \ ; j-^^ife 4|o!;- |al r / 



where (•)^ and (•)^ denote the averaging over the initial amplitudes and initial phases (which can be done independently). 
Our next step will be to calculate the above terms by substituting into them the values of a^^^ and a'^^^ from © and (|10|l 
respectively. 



4.3 Statistical averaging and graphs. 

Let us consider the initial fields 0^(0) = a^"' are essentially RPA as defined above. We will perform averaging over the 
statistics of the initial fields in order to obtain an evolution equations, first for Z and then for the multi-mode PDF. The 
ultimate goal of this exercise is to prove that the wavefield remains of the essentially RPA type over the nonlinear time. 

Let us introduce a graphical classification of the above terms which will allow us to simplify the statistical averaging 
and to understand which terms are dominant. We will only consider here contributions from Ji and J2 which will allow 
us to understand the basic method. Calculation of the rest of the terms, J3, J4 and J5, follows the same principles and 
can be found in Appendix 2. First, The linear in e terms are represented by Ji which, upon using lO, becomes 

Jl = ( n E + ^) {vinamanAi,JL+n + 2t^™a„,a„ Af„<5f+„) d, \ . (24) 

Hereafter we omit, for brevity of notation, the super-script (0) because no other super-scripts will appear from now on. 

Let us introduce some graphical notations for a simple classification of different contributions to this and to other 
(more lengthy) formulae that will follow. Combination V^„5m-|_„ will be marked by a vertex joining three lines with 
in-coming j and out-coming m and n directions. Complex conjugate K^„i5m-i-n '^i^l be drawn by the same vertex but with 
the opposite in-coming and out-coming directions. Presence of aj and dj will be indicated by dashed lines pointing away 
and toward the vertex respectively. ^ Thus, the two terms in formula (|24|l can be schematically represented as follows, 



^This technique provides a useful classification method but not a complete mathematical description of the terms involved. 
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m m 



Let us average over all the independent phase factors in the set {V'}- Such averaging takes into account the statistical 
independence and uniform distribution of variables -0. In particular, (ijj) = 0, {ipi''Pm) = and {'ipi''Pm) = ^i"'- Further, the 
products that involve odd number of tp's are always zero, and among the even products only those can survive that have 
equal numbers of ij/s and ij/s. These ij/s and ij/s must cancel each other which is possible if their indices are matched 
in a pairwise way similarly to the Wick's theorem. The difference with the standard Wick, however, is that there exists 
possibility of not only internal (with respect to the sum) matchings but also external ones with ^'s in the pre- factor H'ipi''' . 

Obviously, non-zero contributions can only arise for terms in which all ij/s cancel out either via internal mutual 
couplings within the sum or via their external couplings to the ^'s in the ^-product. The internal couplings will indicate 
by joining the dashed lines into loops whereas the external matching will be shown as a dashed line pinned by a blob at 
the end. The number of blobs in a particular graph will be called the valence of this graph. 

Note that there will be no contribution from the internal couplings between the incoming and the out-coming lines of 
the same vertex because, due to the 5-symbol, one of the wavenumbers is in this case, which means ^ that V = 0. For 
Ji we have 

Ji = {Ci)4> + (C2)v) 

with 

n 

s 



(Ci)v = > ^ •j + >---<---• 2m 

m 



and 



n 



{C2)^= •j + ntC ^ ^--m2n 

/ 

m 



which correspond to the following expressions, 

^In the present paper we consider only spatially homogeneous wave turbulenee fields. In spatially homogeneous fields, due to momentum 
conservation, there is no coupling to the zero mode fc = because such coupling would violate momentum conservation. Therefore if one of 
the arguments of the interaction matrix element V is equal to zero, the matrix element is identically zero. That is to say that for any spatially 
homogeneous wave turbulence system V^^^ = ^/cl=o/c2 ~ ^/ci/c2=o ~ ^' 
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and 



+ Y.(^2,n + ^)V^ZAlA2mAlZS{^,„, + m^^2r.^-l) J] (25) 

(C2)^= = 2 ^ (A, + ^)i;^-A™A„A, A- 5-+„5(/.,„ + 1)5(m„ - 1)%, ^ n ^(a*') 

+ 2^(A„ + ^)Fl"A2„A^A^"„J(^2„ + l)5(Mn-2) n (26) 

Ti " l^n,2n 

Because of the ^-symbols involving /i's, it takes very special combinations of the arguments ^ in Z{fi} for the terms in 
the above expressions to be non-zero. For example, a particular term in the first sum of H25|l may be non-zero if two /i's 
in the set {fi} are equal to 1 whereas the rest of them are 0. But in this case there is only one other term in this sum 
(corresponding to the exchange of values of n and j) that may be non-zero too. In fact, only utmost two terms in the 
both (|25|) and 12()l) can be non-zero simultaneously. In the other words, each external pinning of the dashed line removes 
summation in one index and, since all the indices are pinned in the above diagrams, we are left with no summation at all 
in Ji i.e. the number of terms in Ji is 0(1) with respect to large N. We will see later that the dominant contributions 
have 0{N'^) terms. Although these terms come in the order, they will be much greater that the terms because the 
limit N oo must always be taken before e — *■ 0. 

Let us consider the first of the e^-terms, J2. Substituting Q into if^ . we have 



>^iVinama,ALJL+n + 2^^" a™a„ A^";,<5;-';„) (y^.a.a, A^,,(5^,+, + 2V;ia^a,A^J^^,))^, 
= {Bi+B2 + B2 + B3)^, (27) 



where 



m K m K m K 

\ J / \ J / \ J / 

Bi= >•-«•••; B2= >--*--s and B3 = v (28) 

/ \ / \ / \ 

n V n v n V 



Here the graphical notation for the interaction coefficients V and the amplitude a is the same as introduced in the 
previous section and the dotted line with index j indicates that there is a summation over j but there is no amplitude aj 
in the corresponding expression. 
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Let us now perform the phase averaging which corresponds to the internal and external couplings of the dashed lines. 
For (-Bi)^ we have 



where 



m K m K 

• ^ • 1/ ^ m 

• ^ • ^ ^ m 



m K 

/ \ 

n V 



m 



+ 2 <...^. ..A , (29) 



/ 



n 



1 iJ?- . - . . ^ . ■ , 



m 

\ 2 ^ 



n 



,2 



J * mn m+n /ti-' 



;A^^„A2-A„A„>l2<5(^„ + l)J(;x„ + l)J(/x,-2) [] 



m 
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We have not written out the third term in H29|) because it is just a complex conjugate of the second one. Observe that 
all the diagrams in the first line of H29|) are 0(1) with respect to large N because all of the summations are lost due to 
the external couplings (compare with the previous section). On the other hand, the diagram in the second line contains 
two purely- internal couplings and is therefore 0{N^). This is because the number of indices over which the summation 
survives is equal to the number of purely internal couplings. Thus, the zero-valent graphs are dominant and we can write 

{Bi)^ = Xl5{f^i) J2 + \-A'j)\ViJ'\AU'Sl+^AlAl[l + 0{1/N^)] (30) 
For (i?2)i/. we have 




n 



> ■ \ = ^ (Aj + XjA"^ — ^-^)V^j„V^'J;A:j,„jA^^(5^^„(5j_^j,Am74„A„yl^ 



m ~m 

1 . / \ ( A - -I- A A ^ T/" A-' 
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n 



/A 

jr \ j=im=i7i=iK j 



m 



2^2 



K 

% m 



V ■ I • 

^ 2m 



2 2 mI 2 



m 



m 
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A 2m A K CK A A 

"mm^2m!y"2m+i/^m^K^i/ 

x,5(^„ + 2)^(^^-l)(5(/i^ + l) J]^ (5(//i) 



/ ,,2 



3rn •-►'■■►■< • = V'(A2m + A2„A2,„ - ^^^f^)V;^™V^2^" 



2m ^ _ 
m 



m 



y m \ 



mm 



:A'Zn^lZmAlA,ra5{^im+m^^3m-l) J] ^(^') 



-m 



^ ' ^ 2 

' ^ "^i — J_ \2 a2 _ JHm^\rr2mrrm A 

\ ^ ^/^I'^Sm ^ '\2m^2m 2y^2 i *^mm ^2m-m ^mm 



2m ^.^^ ,„ ^^2m 

m 



l^m^ — m 

The second term in contains one summation because its graph has one purely internal coupling. This term is N 
times smaller than the largest terms in (-Bi)^ (which have 2 surviving summation indices). All the other terms in 
contain no summation at all because all their dashed lines are coupled externally. 

Similarly, the leading contribution to (B^)^, will be given by the zero-valent graph with the maximum possible number 
of internal couplings (which is equal to 2 in this case). Because of the (5's, there are no graphs with just one internal 
coupling, but there are graphs with all the dashed lines coupled externally. Thus, 
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(B: 



3/V 



m 
j 



y [i + o(i/iv^)] 



2n^(/^') E {^s+>^'Aj)\Vj:::\'\A]l\'SJ\,,,AlAl [1 + 0{1/N% 



Summarising the results of this section we can write for J2: 



m 1 2 I Am 1 2 srm 
jn I I jn I j+n 



AlAl [l + 0(l/iV)]. 



(32) 



(33) 



Thus, we considered in detail the different terms involved in J2 and we found that the dominant contributions come 
from the zero-valent graphs because the have more summation indices involved. This turns out to be the general rule 
that allows one to simplify calculation by discarding a significant number of graphs with non-zero valence. After this 
observation finding the rest of the terms, J3 to J5, becomes a routine task and we therefore move it to the Appendix 2. 



4.4 Equation for Z 

Now we can observe that all contributions to the evolution of Z (namely Ji — J5 , see the previous section and Appendix 
2) contain factor Yli ^(a*j) which means that the phase factors {ip} remain a set of statistically independent (of each each 
other and of A's) variables uniformly distributed on S^. This is true with accuracy O(e^) (assuming that the TV-limit is 
taken first, i.e. 1/A'' <C e^) and this proves persistence of the first of the "essential RPA" properties. Similar result for 
a special class of three-wave systems arising in the solid state physics was previously obtained by Brout and Prigogine 
This result is interesting because it has been obtained without any assumptions on the statistics of the amplitudes 
{A} and, therefore, it is valid beyond the RPA approach. It may appear useful in future for study of fields with random 
phases but correlated amplitudes. 

Let us now derive an evolution equation for the generating functional. Using our results for Ji — J5 in p8(l and (|17|l 
we have 



d 



d^zjo) 

OXmdXn 



+4e^ E A, 



dZjO) 



d^zjo) 

dXjdXndXk 



cc. 



(34) 



Here partial derivatives with respect to A/ appeared because of the Ai factors. This expression is valid up to 0(e'*) 
and 0{e^ /N) corrections. Note that we still have not used any assumption about the statistics of A's. This is a linear 
equation: as usual in statistics we traded nonlinearity for higher dimensions. The last term here "spoils" the separation 
of variables and, therefore, puts a question mark on the independence of variables {A} from each other on the nonlinear 
time. 

Let us now N 00 limit followed by T '--^ 1/e ^ cxd (we re- iterate that this order of the limits is essential). Taking 
into account that lim E{0,x) = T{tt6{x) + iPi^)), and lim \A{x)\^ = 2ttT6{x) and, replacing (Z(T) - Z(0))/r by Z 
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we have 



-2Xi 



_5 5_ 

5^Z 



sz 



SXjSXnSXr 



■ J' dkj dkuidkfi 



(35) 



Here variational derivatives appeared instead of partial derivatives because of the TV — > oo limit. 



4.5 Equation for the PDF 

Taking the inverse Laplace transform of (|35|l we have the following equation for the PDF, 



V = - 



(36) 



where Fj is a flux of probability in the space of the amplitude s, , 



-2V 



-2(|^^"„|2<5(c.^„)<5;.V„ - 2\VlS5{^Ln)5L+n)j-{-W^V)] dk^dK 



(37) 



This expression can be simplified to 



F, 



{{\V,UH{u;USL+n + 2\V^lfS{u;^JS^^JsnSr 



+2Vi\Vpj^6{co^l^)dJ+„, - |K^,„P<5(a>^„)<5^+„)s„ 
+2{\Vpj^d{u^Jd^+,^ - 2\ViJ'S{Loi,JSl+Js„. 



6V 



5V 



J' dkffidJxiji 



(38) 



This equation is identical to the one originally obtained by Peierls ^H] and later rediscovered by Brout and Prigogine [Tl)] 
in the context of the physics of anharmonic crystals. Zaslavski and Sagdeev jl7j were the first to study this equation in 
the WT context. However, the analysis of JX was restricted to the interaction Hamiltonians of the "potential 

energy" type, i.e. the ones that involve only the coordinates but not the momenta. This restriction leaves aside a great 
many important WT systems, e.g. the capillary, Rossby, internal and MHD waves. Our result above indicates that the 
Peierls equation is also valid in the most general case of 3-wave systems. 

Here we should again emphasise importance of the taken order of limits, — s- oo first and e — > second. Physically 
this means that the frequency resonance is broad enough to cover great many modes. Some authors, e.g. [1511161 [T7| . leave 
the sum notation in the PDF equation even after the e ^ limit taken giving 5(ijj",^). One has to be careful interpreting 
such formulae because formally the RHS is nill in most of the cases because there may be no exact resonances between 
the discrete k modes (as it is the case, e.g. for the capillary waves). In real finite-size physical systems, this condition 
means that the wave amplitudes, although small, should not be too small so that the frequency broadening is sufficient to 
allow the resonant interactions. Our functional integral notation is meant to indicate that the N oo limit has already 
been taken. 
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5 Approximate independence of the amplitudes. 



Variables Sj do not separate in the above equation for the PDF. Indeed, substituting 

'p{N,a) _ p(a) p{a) p(a) ^^g-j 

into the discrete version of H38I) we see that it turns into zero on the thermodynamic solution with pj"'' = LOj cxp{—LUjSj). 
However, it is not zero for the one-mode PDF pj"'' corresponding to the cascade-type Kolmogorov-Zakharov (KZ) spec- 
trum Tij^, i.e. pj-"-* = (l/rij^) eKp{—Sj /rij^) (see next section), nor it is likely to be zero for any other PDF of form H39() . 
This means that, even initially independent, the amplitudes will correlate with each other at the nonlinear time. Does 
this mean that the existing WT theory, and in particular the kinetic equation, is invalid? 

To answer to this question let us differentiate the discrete version of the equation (|35|l with respect to A's to get 
equations for the amplitude moments. We can easily see that 

d,{{AlAl)~-{Al)(Al))=Oie^) (ji,j2 e S^) (40) 

if (Aj^Aj^Aj,^) = (A^J (A^J (A^J (with the same accuracy) at t = 0. Similarly, in terms of PDF's 

{pj':Si^n,^3.) - Pj:\^3.)Pj:\^3.)) = Oie') (ji, J2 e Bn) (41) 

if ^£.3..J^.- ^.3, = P^:\sn)P^:\^n)P^:\^3s)Pi:\s3.) at * = O. Here Pi!:S,,,,M,, s,,, s,,, s,J, P^^^(.,„ 
and Pj'^\sj) are the four-mode, two-mode and one-mode PDF's obtained from V by integrating out all but 3,2 or 1 ar- 
guments respectively. One can see that, with a accuracy, the Fourier modes will remain independent of each other in 
any pair over the nonlinear time if they were independent in every triplet at t = 0. 

Similarly, one can show that the modes will remain independent over the nonlinear time in any subset of M < N 
modes with accuracy M/N (and e^) if they were initially independent in every subset of size M + 1. Namely 

Pj'!^:!...,M i^n , , ^.M ) - P^:^ i^'^n )P\:^ i^. )■■■ (^.m ) = 0{M/N) + 0{e^) 

(ji,j2,...,jAf e 6jv) (42) 

:f p(M+l,a) _ p(o) p(a) p(a) „j. ^ _ n 

Mismatch 0{M/N) arises from some terms in the ZS equation with coinciding indices j. For M = 2 there is only 
one such term in the A^-sum and, therefore, the corresponding error is 0{1/N) which is much less than O(e^) (due to the 
order of the limits in N and e). However, the number of such terms grows as M and the error accumulates to 0{M/N) 
which can greatly exceed O(e^) for sufficiently large M . 

We see that the accuracy with which the modes remain independent in a subset is worse for larger subsets and that 
the independence property is completely lost for subsets approaching in size the entire set, M ~ A^. One should not 
worry too much about this loss because N is the biggest parameter in the problem (size of the box) and the modes will 
be independent in all Af-subsets no matter how large. Thus, the statistical objects involving any finite number of modes 
are factorisable as products of the one-mode objects and, therefore, the WT theory reduces to considering the one-mode 
objects. This results explains why we re-defined RPA in its relaxed "essential RPA" form. Indeed, in this form RPA 
is sufficient for the WT closure and, on the other hand, it remains valid over the nonlinear time. In particular, only 
property (|40|l is needed, as far as the amplitude statistics is concerned, for deriving the 3-wave kinetic equation, and this 
fact validates this equation and all of its solutions, including the KZ spectrum which plays an important role in WT. 

The situation were modes can be considered as independent when taken in relatively small sets but should be treated 
as dependent in the context of much larger sets is not so unusual in physics. Consider for example a distribution of 
electrons and ions in plasma. The full A^-particle distribution function in this case satisfies the Liouville equation which 
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is, in general, not a separable equation. In other words, the A^-particle distribution function cannot be written as a 
product of N one-particle distribution functions. However, an M-particle distribution can indeed be represented as a 
product of M one-particle distributions if M ^ Njj where Njj is the number of particles in the Debye sphere. We see 
an interesting transition from a an individual to collective behaviour when the number of particles approaches Njj . In 
the special case of the one-particle function we have here the famous mean-field Vlasov equation which is valid up to 
0{1/N£,) corrections (representing particle collisions). 



6 One-mode statistics 

We have established above that the one-point statistics is at the heart of the WT theory. All one-point statistical objects 
can be derived from the one-point amplitude generating function, 

which can be obtained from the A''-point Z by taking all /^'s and all A's, except for Aj, equal to zero. Substituting such 
values to (|35(l we get the following equation for Za, 

^ = A,r;,Z. + (A,^ry,-A,7.)^, (43) 

where, 

77, - 4^6^ j {\VL?5iJ{ujlJ + 2\V;r?5]tS{0JTi)) ^Mn^ dkidk^, (44) 

7, = 87re2 j (\Vil\HlJ{ujlJn„, + \V^\Hfi5{uJ){ni - n„S) dhdk^. (45) 

Correspondingly, for the one mode PDF Pa{sj) we have 

dPa dF , , 

F = -s,{lPa^-^,^). (47) 

Equations H43() and (|46|l where previously obtained and studied in [201 in for the four-wave systems. The only difference 
for the four-wave case was different expressions for t] and 7. For the three-wave case, equation for the PDF was not 
considered before, but equations for its moments were derived and solved in |19| . In particular, equation for the first 
moment is nothing but the familiar kinetic equation h = —771 + r\ which gives r] = 771 for any steady state. This, in turn 
means that in the steady state with = we have ^j"^ — i^/nj) exp(— Sj/?T,j) where rij can be any steady state solution 
of th kinetic equation including the KZ spectrum which plays the central role in WT [Sm. However, it was shown in 
pn| that there also exist solutions with F 7^ which describe WT intermittency. 



with is a probability flux in the s-space, 



7 Phase statistics. 

Importantly, RPA formulation involves independent phase factors i/j = e*^ and not phases (p themselves. Firstly, the 
phases would not be convenient because, as we will see later, the mean value of the phases is evolving and one could 
not say that they are "distributed uniformly from — tt to vr". In fact, we will also see that the mean fluctuation of the 
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phase distribution is also growing and they quickly spread beyond their initial 27r-"wide interval. But perhaps even more 
important, 0's build mutual correlations on the nonlinear time whereas Tp^s remain independent. This will be shown later 
in this section, but we would like first to give a simple example illustrating how this property is possible due to the fact 
that correspondence between </) and is not a bijection. 

Let TV be a random integer and let ri and r2 be two independent (of N and of each other) random numbers with 
uniform distribution between — tt and tt. Let 

4>i^2 — "^ttN + ri^2- 

Then 

(<^i,2) =27r(Ar), 

and 
Thus, 

- {<Pi){h) = 4^'((A^') - W) > 0, 
which means that variables and 02 are correlated. On the other hand, if we introduce 

then 

(^1.2) = 0, 

and 

= 0. 

(V'l^2)-(V'l)(^2)=0, 

which means that variables ^1 and ^2 are statistically independent. In this illustrative example it is clear that the 
difference in statistical properties between </) and arises from the fact that function "0(0) does not have inverse and, 
consequently, the information about N contained in </> is lost in tp. 

This illustration, although simple, captures the property that actually happens in reality as we will show below. Let 
us use the following expression for the phase 

0j = 5 In aj. 

Substituting ||HJ| and Taylor-expanding of logarithm in e one gets 



(1) _ 



= 31n(af + eaf^ + e'af) = cl>^"^ + ecf,^'^ + £^0(2) (48) 
where 

(49) 
(50) 

(51) 

Now let us perform averaging over the statistics of factors -0'"^ . As usual, the surviving terms are those in which all -0^"^ 's 
cancel out due to their pairwise matchings. This is possible only if the number of '0(''''s is equal to the number of '0(''''s 
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in the products defining these terms. Easy to see that the e term involves three -^^'^^'s and therefore its average is zero. 
Therefore, 

(MT)) - (0f > = (-^ {4^^'y+4^^') (52) 

Let us consider 



\m.n.K,y 



Here, there are two terms with equal number of i/i^^^'s and ^/j^'^^'s but all couplings of index j to any other index give zero 
because V = if one of its wavenumbers is zero. Thus, {iP'^po!'^^)'^'^ = 0- The other term, (yO'f^o,f^'^ has already been 
calculated before when evaluating J3. We have 



(2)-(0)^ 

a) 'a) ' 
1 3 / ^ 



4 



Let us take limits N ^ 00 and T ^ 00 and replace {4>{T))^ — (0(O))^/T by (^)^. We get 

= ojnl, (53) 

where conl is the nonlinear frequency correction given by 



= 4e' / \\VL\'-P (^) - IVJ^l'-P (^] ST+ni^l - Al 



A^- dkmdkn (54) 



Here P{x) denotes the principal value of the integral. Averaging over the amplitudes, we have 

{(f>j) = {ojnl), 

where {lonl) is the amplitude-averaged nonlinear frequency correction 



{uj^^) = 4e^ / {^) 6L+n - \Vj::\'P l;^] SJ!,^{nm - n„) 



We can see that the mean value of the phase is steadily changing over the nonlinear time and, therefore, it would be 
incorrect to assume that the phase "remains uniformly distributed from — tt to tt" even though this could be true for 
t = 0. This is one of the reasons why we formulate RPA in terms of and not 4>. Indeed, i/j was shown above to stay 
uniformly distributed on the unit circle over the nonlinear time. 

The other reason is that, strictly speaking, ^'s do not stay de-correlated where as ip's do (as shown before). We 
already saw in the beginning of this section that this situation is possible due to the fact that the map 4) ^ ijj = e"^ 
is not a bijection. Let us now study such a buildup in statistical dependence of the phases, let us consider correlator 
^3,k = {{4>j - {(f>j)){^k - {(t>k))) = {4>j(t>k) - 4>k)- At time T we have 

^3,k{T) = + e^^) + £^4^2 (56) 
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where 



•r-(O) 


= (<^;% 


4"^> - ( 


4"' 








= 










^(2) 


= (0^. 











(57) 



Here, we have taken into account that, as we showed earher, {(j^^p) = 0. Let us consider the e-term ^j^, e.g. 



r(i) 



(l)-(0)^ ,(0) 



fc,m,n 



(58) 



In this expression, we have a factor cpk which enters directly and not via the combination tpi- = e^^'' . Potentially, this 

could greatly complicate the situation because to objects like {ijjk4'k) knowledge of the statistics of ■;/' is not sufficient and 
one needs the full PDF of (pk- Fortunately, however, this does not cause problems here because, no matter what index is 
matched to k, matching of the two remaining indices results in F = 0. Therefore, the contribution of the e-terms is nill. 

(2) 

Let us now consider the .Fj starting with 




(2) 

i 

(0) 

j 



2A] 



(59) 



We see that the square bracket on the RHS involves an even number (four or six) of ijj's in each term. Thus, in 
order for these terms to survive these ijj's must cancel out which is possible when their indices match in a pairwise 
way. But this means that index k (of 4>k) docs not match to any of the indices of i/j's and, therefore, the averaging of 
(j>k can be taken separately because it is statistically independent of all other phase factors.^ Thus, we conclude that 



a(2)w^(o)\ 



= {<t>'f')i!{4''k')4> and these terms drop out of J-'^^l- The remaining term in T', 



(2) 
j-k • 



,(2) 

j,k 



+ 




^\nu-k/\^' 

si 



k{j-k)\ ^j-k 



3{k-i)\ \^3{k-i)\ ^k-j+'^\^jk 



+ C.C 



l,m ^ 



1 l2^J , 9|T/m|2| A m|2 cm 



+fe|2| Ai+fe|2 
^jk 



"^j+k 



(60) 



We can now average over the amplitudes and take limits N ^ 00 and e — > and write 



,j+k\ 



+ 



j \\yL?W5i+m + 2\vjr\H{ujj})sr+, 



niTim dkidk 



(61) 



Presence of the 1-st term on the RHS indicates that the phases of the j-th and the fc-th modes get correlated on the 
nonlinear time. This correlation is week in a sense that J^j^k has a sharp peak at j = k but the integrated contribution 
of all j ^ k'ls of the same order as the value at the contribution of the j = k peak and, therefore, could cause a problem 
should one tried to build RPA based on the statistics of ^'s rather than i/)'s (which remain de-correlated). 

^ There is of eourse also a possibility that k couples simultaneously to both indices in a pair, but this contribution contains ~ 1 /N less terms 
and, therefore, should be ignored. 
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Let us consider a special case of ()m|l for j — k which is interesting because it allows one to calculate the dispersion in 
phases, 

We have 

ijk^Vk/nk, (62) 

where rjk is defined in (|44(l and Uk = (la^p). One can see that the RHS here is always positive and, therefore, the phase 
fluctuations experience an unlimited growth. On stationary spectra, this growth is ^ \/t which corresponds to cr ~ i. 
Recall that the mean value of the phase is also changing in time with the rate lon l and on stationary spectra this change 
is linear in time. 



8 Discussion 

In the present paper, we considered evolution of the full N-mode objects such as the generating functional and the 
probability density function for all the wave amplitudes and their phase factors. We proved that the phase factors, being 
statistically independent and uniform on initially, remain so over the nonlinear evolution time in the leading order in 
small nonlinearity. If in addition the initial amplitudes are independent too, then they remain so over the nonlinear time 
in a weak sense. Namely, all joint PDF's for the number of modes M <^ N split into products of the one-mode densities 
with 0{M/N) and O(e^) accuracy. Thus, the full iV-mode PDF does not factorise as a product of N one-mode densities 
and the Fourier modes in the set considered as a whole are not independent. However, the wave turbulence closure only 
deals with the joint objects of the finite size M of variables while taking N ^ oo limit. These objects do factorise into 
products and, for the WT purposes, the Fourier modes can be interpreted as statistically independent. In particular, the 
derivation of the kinetic equation for the energy spectrum deals only with the 1-mode and the 2-mode distributions and 
is, therefore, justified by the results of the present paper. Generally speaking, our results reduce the leading-order WT 
problem to the study of the one-mode amplitude PDF's and they validate the generalised RPA technique introduced in 
[T^ 1201 ■ Such a study of the one- mode PDF and the high-order momenta of the wave amplitudes was done in |2n| . 
It was shown, in particular, that anomalous probabilities of large wave amplitudes can appear in the form of a finite-flux 
solution in the amplitude space caused by a wave-breaking amplitude cutoff. The reader is referred to these papers for 
the discussion of the WT intermittency. 

Although our results indicate that correlations between 2 or more (but ^ N) modes do not appear in the leading 
(i.e. e^) order for the three-wave systems, they definitely appear as corrections in the next (i.e. e^) order. Our paper 
is concerned with the main order statistics only in which the main evolution happens in the 1-mode objects, e.g. the 
1-mode amplitude distributions. For study of the multi-mode correlations developing in WT in the next order in e the 
reader is referred to papers |21[ [221 ■ 

We have also considered correlators of the phase and we showed the relation between the statistical properties of the 
phase (j) and the phase factors ip = e**^. We showed that the mean of (j) and its fiuctuations about the mean grow in time 
and, therefore, there exist no 27r-wide interval in which the phase would remain uniformly distributed. Moreover, phases 
(j) become correlated at different wavenumbers that lie on the resonant manifold. These properties make the phase cj) an 
inappropriate variable for formulating the RPA method of WT description. On the other hand, our work shows that the 
phase factors = e*"^ do remain statistically independent and uniform on which makes them the right choice for the 
RPA formulation. 

The present paper deals with the three-wave systems only. The four-wave resonant interactions are slightly more 
complicated in that the nonlinear frequency shift occurs at a lower order in nonlinearity parameter than the nonlinear 
evolution of the wave amplitudes. To build a consistent description of the amplitude moments one has to perform a 
renormalisation of the perturbation series taking into account the nonlinear frequency shift. This derivation will be 
published separately, whereas here we just announce its main result, the 4-wave generalisation of the Peierls equation for 
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the PDF. It has the same continuity equation form (|36|1 but now the probabiHty flux is 



123 J 1 z a 

where W^^ is the 4-wave ineraction coefficient and o)^" — uj'i^^ + il/ + fi^ — — 17,^ with fi; = 2e^^ VF/^ti^ being 
the nonhnear frequency shift. As wee see this equation is even more compact than its 3-wave analog. In addition 
to the derivation of this equation, we will also analyse its properties and consequences for the mode correlations and 
intermittency in 4-wave turbulent systems. 



9 Appendix 1 



Let us obtain Z{T) in terms of the series in small nonlinearity up to the second order in e. As an intermediate step, we 
first consider separately the amplitude and the phase ingredients of Z and substitute the e-expansion of a from (jS)) into 
their expressions. 



O J I J ' J J ' 



A,|a("'|2+eA,(a(l'a<.°)W.l'a<."')+£2Aj|a(l'P + (a<.=)a(°'+a^ 



A,A 



/ 2 \ 

J (1 + eaij + e a2j), 



(64) 



and 



(0)^3 



^f^^^(l + 6/?i,+e2/?2,). 



(0)Mj 



^ 3 



2 2 V 2 



/ (2) 

I ay a 



(0) 



1 ^2 S 

2 af ) 2 ,(0) 



^3 

,(0) 



,(0) 



Mi 
4 



^ _ 1 



(1) \ 2 



2 2 V 2 



1 

my 



(1) - 
W 



2|„(1)|2 



(0)2 



(6 



where aij,a2j,l3ij and j32j denote the linear and quadratic contributions into the amplitude and phase parts of Z 
respectively. 



ay 

P2j 



(A, + A|Af >f )p + A,(af '4" ) + af )af ) + |(af ^af V + 



2yl 



(0)2' 



3 



2A 



(0)2 



(1)^ 



2|„(1) 12 



4A 



(0)2 



(66) 
(67) 
(68) 

(69) 
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Substituting expansions ()64|l and (|65f) into the expression for Z, we have 



Z{T) 



i2n) 



N 



n 



ai 



1 /TTr^-4i°'%,.(0)M, 



j<fc 



For parts /i,/2 and /a in the above expression we have, 



E 



2 ^ 4A! V 2 V 



2A2 



2 ^ 4A^{2 ^ J 



2A2 



Exploiting the property Z{X, — /i} = Z{A,/i} we can write 

Z{\,fi}=X{X,fi} + X{X,~fi}. (70) 

At t = T we have for X{X, ^} 

X(r) - X(0) + (27r)2^ / n e^'l^'°'l'[eJi+e2(J2 + J3 + J4 + J5)]\ , (71) 

\\\1\\<N I A 



where 



2\a 



2|a 



(0)|2- 
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J4 



+ 



2 ■ 4|af |4 



+ 



2|o' 



(0)|2 



/ (1)_(0)n2 

(a) ^a} 



(74) 
(75) 



= I (n W"^''' E A. A.(a«af + a^.^^af )ai^)ar + (A, 

\ I jy^k 



Mfc ^J1):,(0) -(l)„(0)x (l)-(O) 



4|a 



where (•)^ and (•)^ denote the averaging over the initial amphtudes and initial phases respectively. We remind that such 
individual averages are possible because the amplitudes and the phases are statistically independent from each other at 
t = 0. 

10 Appendix 2 
10.1 Calculation of J3 



^3 = nv^^"'"E(A.+^)«. 



2^2 



^ [2y^„ (-V;™a„a«a,ii;[< 



(77) 



The terms to be averaged here can be drawn as 



+ 



V n 



K 3 



+ 



+ ;• • •►• < 

/ m ^ 



V n 



(78) 
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Let us now average over the random phases. Again, the leading order terms will be given by the diagrams with the 

largest number of internal couplings. They will arise from the VV terms (the 2nd, 3rd and the 6th graphs) because they 
allow 2 internal couplings in each of them. There are also possibilities to have one internal and two external couplings of 
the dashed lines, - such terms will give a 0{1/N) correction to the leading order. Therefore, 
j 3 3 




[l + 0(l/iV)] 



n n m 



x[l + 0(l/iV)] 



(79) 



10.2 Calculation of J4 
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Graphically, the 4 terms to be averaged in this expression are 



(80) 
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Note that there is no dotted lines in these graphs because for each summation index there is a corresponding wave 
amplitude present. As a consequence, the rule for the number of surviving summations is somewhat different from what 
we had so far. Namely, the number of the summation indices after the phase averaging is one less than the number of the 
purely internal couplings. Easy to see that the phase averaging of the above terms always leads to an external coupling 
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of the dashed lines j which removes the j summation. Moreover, no more than one purely internal coupling of the dashed 
lines is possible in any of these graphs. ^ Thus, J4 contains no summation at all and is only a 0{1/N'^) correction to the 
main terms in J2 and J3. 

10.3 Calculation of J5 

Expression for J5 seemingly involves a great number terms. However, this number can be dramatically reduced by the 
following speculation. Just as in J4 there is no dotted lines in the graphs involved in J5 because for each summation index 
there is a corresponding wave amplitude present. Thus, we have the same rule for the number of summations surviving 
the phase averaging (i.e. one less than the number of internal couplings) . In order to be of the same order as the leading 
terms in J2 and ./-j. wo must have 3 purely internal couplings and, therefore, no external couplings. This is only possible 
when the number of dotted lines directed to the vertices is equal to the number of them pointing away which is true for 
the VV terms but not true for the VV and VV terms. Thus we will only consider the VV terms. Further, the fact that 
there is no external couplings means that such terms are only non-zero when all /i's are zero. Thus, there will be no 
contribution from the second part of J5 which has a Hk pre-factor. 

where the VV terms to be averaged here are 
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m ^ K 
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®The only possibility of the double internal coupling would be in the j^t graph via joining m with v and K with n, but this would mean 
j = because of the 5-symbols and, therefore, this term is nill. 
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By coupling the dashed Hnes we have in the leading order 



(Ci)^ = 



k 
n 



2 I 

\ 



j 
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Term C4 does not survive the averaging because j ^ fc and a triple internal coupling is not possible. Summarising, 



J5 = 2llH^^l) E ^^■^'^ 



AjAiAi [l + 0(l/iV)]. 



(82) 
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